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ABSTRACT 

A perturbative analysis is used to investigate the effect of rotation on the 
instability of a steady accretion shock (SASI) in a simple toy-model, in view of 
better understanding supernova explosions in which the collapsing core contains 
angular momentum. A cylindrical geometry is chosen for the sake of simplicity. 
Even when the centrifugal force is very small, rotation can have a strong effect on 
the non-axisymmetric modes of SASI by increasing the growth rate of the spiral 
modes rotating in the same direction as the steady flow. Counter-rotating spiral 
modes are significantly damped, while axisymmetric modes are hardly affected 
by rotation. The growth rates of spiral modes have a nearly linear dependence on 
the specific angular momentum of the flow. The fundamental one-armed spiral 
mode (m = 1) is favoured for small rotation rates, whereas stronger rotation rates 
favour the mode m = 2. A WKB analysis of higher harmonics indicates that the 
efficiency of the advective-acoustic cycles associated to spiral modes is strongly 
affected by rotation in the same manner as low frequency modes, whereas the 
purely acoustic cycles are stable. These results suggest that the linear phase of 
SASI in rotating core-collapse supernovae naturally selects a spiral mode rotating 
in the same direction of the flow, as observed in the 3D numerical simulations 
of Blondin & Mezzacappa (2007). This emphasizes the need for a 3D approach 
of rotating core-collapse, before conclusions on the explosion mechanisms and 
pulsar kicks can be drawn. 



Subject headings: hydrodynamics — shock waves — instabilities — supernovae: 
general 
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Introduction 



Despite extensive studies, the explosion mechanism of core-collapse supernovae is still 
elusive. According to the delayed explosion scenario, the shock is first stalled at a distance of 
a few hundred kilometers, and then revived after neutrinos diffuse out of the proto neutron 
star. Unfortunately, numerical simulatio ns suggest that neutrino h eating may not be efficient 
enough, at least in spherical symmetry (ILiebendorfer et al. 1120051 ). 



Recent studies have shown that the spherical stalled shock is unstable against non ra- 
dial perturbations with a low degree I = 1, 2 even if the flow is conv ectively stable. This 



result was demonstrated using axisymmetric numerical simulations (iBlondin et al.l 12003 



Blondin fc Mezzacappa 



2006 



Scheck et al 



2006 



2008 



Ohnishi et al.l l2006h and linear sta- 



bility analyses (IGalletti &: Foglizzo 1120051 ; iFoglizzo et al.l 120071 ; lYamasaki k, Yamadal 120071 ) . 



Some numerical simulations have shown that this hydrodynamical instability, often called 
SASI, may assist the revival of the shoc k and trigger a success ful explosion, powere d either by 
neutrino heating (jMarek Sz Jankall2007 ) or by acoustic waves ( Burrows et al. hood ). Some ob- 
served properties of young neutron stars may a lso be the consequenc es of SASI, such as their 
distribution of velocities (jScheck et al.l 12004 . 120061 ) or their spin (IBlondin fc Mezzacappa 
20071 : iBlondin fc Shawl l2007h . 



Until now, most studies of SASI have assumed that the unperturbed flow is purely ra- 
dial and not rotat ing. Since the angular momentum of massive stars is likely to be large 
( IHeger et al.l 120051 ). it is desirable to understand how the properties of SASI are affected by 
rotation. In this Letter, the effect of rotation on the linear stage of SASI is investigated 
usi ng a perturbative analys i s in or der to shed light on one of the surprising results observed 
by IBlondin fc Mezzacappa! (120071 ) in their 3D numerical simulations: the development of 
SASI seems to systematically favour a spiral mode rotating in the same direction as the 
accretion flow. As a consequence of momentum conservation, this mode diminishes and may 
even reverse the angular momentum acquired by the proto-neutron star from the station- 
ar y flow. Incidentally, the prese nt linear study does not address another surprising result 
of IBlondin fc Mezzacappal (120071 ). that a spiral mode of SASI always dominate the axisym- 
metric mode even without rotation. Following an approach similar to IFoglizzo et al.l (120071 ) 
(hereafter FGSJ07), we first compute the eigenfrequencies by solving accurately a boundary 
value problem between the shock surface and the accretor surface; in a second step, we use the 
same WKB method as in FGSJ07 to measure the stability of pur ely acoust i c and advective- 
acoustic cycles in this region. This approach is different from lLaming I (120071 ). which is 
based on the approximate derivation of a dispersion relation. Rather than the complexity 
of de scribing the non-spherical shape of a shock deformed by rotation (lYamasaki fc Yamada 
20051 ). we have chosen, as a first step, to solve the much simpler problem of a cylindrical 



- 3- 



accretion shock. This flow is simple enough to allow for a complete coverage of the parame- 
ter space and a physical insight of the main effects of rotation on SASI. Once characterized, 
these effects can be transposed into the more complex geometry of a rotating stellar core. 



Formulation 



Following a similar desire of simplification, iBlondin &: Shawl (120071 ) have limited their 
2D simulation domain (r,<p) to the vicinity of the equatorial region (6 = it/ 2) of a non- 
rotating spherical flow, where the poloidal velocity of a symmetric mode is vg = 0. Their 
calculation neglected the perturbation of density induced by the term (p/r)(dvg/d9) in the 
equation of continuity, implicitly decoupling the structure of perturbations in the vertical 
direction from their structure in the equatorial plane. In order to study the effects of rotation 
in a simple setup and maintain a self consistent set of equations, we have chosen to adopt a 
cylindrical geometry invariant along the z-axis of rotation. The accretor and the shock are 
thus cylindrical in our toy-model. 

Except for the geometry of the flow and the presence of rotation, the assumptions are the 
same as in FGSJ07. In the equations describing the flow (Appendix A), we have assumed 
that i) the free-falling supersonic flow is cold, ii) gravity is Newtonian and self-gravity is 
neglected, iii) the shock is adiabatic, iv) neutrino heating is neglected and the cooling rate 
per volume is approximated as £ oc pP~ a P a with a = 3/2, and (3 = 5/2, where p is the 
density and P the pressure, v) the accreting material is described by a perfect gas with a 
uniform adiabatic index 7 = 4/3, vi) the condition Ai = is imposed at the inner boundary, 
where M. = —v r /c is the Mach number associated to the radial velocity v r and adiabatic 
sound speed c, vii) the flow is inviscid and its specific angular momentum L is conserved. 

We have adopted values of the parameters typical for the core-collapse problem: the ra- 
dius of the inner boundary is r* = 50 [km] and the gravitational potential is $ = — 1.3GM@/r, 
where G is the gravitational constant and M Q the solar mass. By adopting the same gravi- 
tational potential as for a spherical accretor, the Bernoulli equation is unchanged. 

The stability of the accretion flow is investigated for various values of the specific an- 
gular momentum L, measured by the rotation frequency f p = L/(2nr^) extrapolated at a 
radius r p ~ 10 [km] by reference to young pulsars. The range of rotation rates considered 
corresponds to < f p < 10 3 [Hz]. We have chosen to compare instability timescales in cylin- 
drical models with different rotation rates but identical shock radii and mass accretion rates, 
in order to keep geometrical factors constant. This is made possible by adapting the inten- 
sity of cooling accordingly, by a modest amount (< 15%) over the range of rotation rates 
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considered. We observe that the advection time from the shock to the accretor is hardly 
affected by rotation (a few percent). This is because the centrifugal force L 2 /r 3 is much 
smaller than gravity —d&/dr: 



^ ~ 4.6x 10- 2 f^) (-A—) 2 . (1) 



r 3 \d<S>/dr\ " \rJ \10 3 [E 

The perturbations superimposed upon the steady solution are proportional to exp{— i(cut — 
m6 — k z z)}. We adopt the variables 5S, 5q, 5f, 5h defined by 

is s J- «£_«£, (2) 

7 — 1 c p 

Sq = s( [ r —dA, (3) 



pv. 



r 



L 2 

5f = v r 5v r H — 5vg H c5c — Sq, (4) 

r 7 — f 

r , <5r r 5p 

5/i = — + — , 5 

v r p 

where 5 denotes the Eulerian perturbation and v e is the azimuthal velocity. Then the differ- 
ential system describing the perturbations becomes particularly compact: 

dr v r (l — M 2 ) { u c 2 

+ [l + (7-l)X 2 ]^-§}, (6) 

dr v r (l - M 2 ) [c 2 oo c 2 ) 

d5S iu' rri . / C \ , nN 
+ 5 -— , (8) 



dr f r \ Pf 

d5g iu' . ( C , 

+ <J — , (9) 



dr r r \ pr r 

where /x and a;' are defined by 

„ 2 = 1 + (10) 

w = • ( n ) 

These equations are solved by imposing the Rankine-Hugoniot relations for the perturbed 
quantities, which are written as, 

= tv rA A( 1 , 12 

OJ V V r A J 
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Sh sh = -i— AC(1-^J, (13) 

V r ,sh V V r,l J 

5S sh ■ U ' V r,l K/ , ( ^r.shV >C sh -£iAC 



7 C^ h V V r,l J Psh^r,sh C^ h 



^f^^+^-^y (14) 



5g sh = ^ac, (15) 

Psh^V,sh 

where the subscripts 'sh' and '1' refer to the values just below and above the shock, respec- 
tively; A( is the radial displacement of the shock surface. Since we have assumed that the 
flow above the shock is cold, the cooling rate C\ = 0. In addition to these equations, we 
impose the condition that the radial velocity perturbation vanishes at the inner boundary 
(5v r = 0). The derivations of the basic equations and the boundary conditions are shown in 
the Appendix A. 

When rotation is neglected (L = 0, u' = u), we remark the formal resemblance between 
the above formulation and the formulation of FGSJ07 describing a spherical flow. The only 
difference is the expression for the parameter /i 2 in Eq. (FIU|) . where (m 2 + r 2 k 2 z ) replaces 
1(1 + 1), and a geometrical factor 2 in the boundary condition for the entropy perturbation 
(v r shV r l replaces 2v r ^v r ^ in Eq. [Tj 



For a small value of the angular momentum L, we also remark that the effect of the 
centrifugal force L 2 /r 3 on the stationary flow is quadratic, and so is its effect on the entropy 
perturbation in Eq. (|14p . The only first order effect of rotation on the differential system 
satisfied by Sf, Sh, SS, Sq, is the Doppler shift described by uj' in Eq. (TTTT) . 



3. Results 

Given the resemblance of formulations, let us first check that the stability properties of 
the cylindrical flow without rotation resemble those of the spherical flow studied by FGSJ07 
despite the different geometry. In the spherical problem, the axisymmetric mode m = and 
the spiral modes ±m of a given degree / have exactly the same growth rate (FGSJ07). In 
the cylindrical flow, the spiral modes ±m are also degenerate without rotation. A numerical 
resolution of the eigenfrequencies shows that the instability is dominated by the mode m = 
±1 if rsh/r* > 2, and by a larger \m\ for smaller shock radius, exactly as observed in the 
spherical flow (Fig. 6 of FGSJ07). The growth rate of the axisymmetric mode (m = 0, 
k z > 0) however, is expected to differ from the spiral modes (m > 0, k z = 0) in a cylindrical 
flow. For example if r s h/r* = 5, the most unstable axisymmetric mode is twice as slow as 
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the most unstable spiral mode. 

The dependence of the growth rate on the specific angular momentum is shown in 
Fig. Q] for the spiral modes m = ±1,±2 in a flow where r^/V* = 5. The growth rate 
of the modes rotating in the same direction as the flow (m > 0) is increased by rotation, 
whereas the counter-rotating modes (m < 0) are stabilized. The increase of the growth 
rate is almost proportional to the specific angular momentum. In a rotating flow with 
L ~ 200 x 27r • 10 12 [cm 2 /s] (f p = 200 [Hz]), the growth rate of the fundamental mode m = 1 
is twice its value in a non rotating flow of same size. The strong effect of rotation on the 
growth rate of SASI does not seem to depend on the presence of a corotation radius r co , 
defined by Re(u;) — mLjr 2 CQ = 0, also displayed in Fig. [TJ We also investigated the effects of 
rotation on the axisymmetric modes (m = 0) for various values of the vertical wave number 
k z , and found that their growth rates are hardly affected, by less than one percent of |f r ,shlA"sh 
in our study (Fig. |2j). A global overview of the parameter space of the cylindrical SASI is 
displayed in Fig. [31 which indicates the azimuthal wavenumber m of the most unstable mode 
for a wide range of shock radii and specific angular momentum. The asymmetric one-armed 
spiral mode is unstable in most of the parameter space, and always more unstable than 
without rotation. The growth rate of the mode m = 2 exceeds that of the mode m = 1 (by 
less than 10% in the example of Fig. 1), as the specific angular momentum is increased. 



4. Discussion 



4.1. Instability mechanism 



As underlined in Sect. 3, the dynamical effect of the centrifugal force on the stationary 
flow is modest. We anticipated in Sect. 2 that the only linear effect of angular momentum 
is a Doppler shift of the eigenfrequency u' = uo — mQ(r), where Q(r) is the local rotation 
frequency. This leaves the axisymmetric mode m = unaffected and explains the relative 
insensitivity of its growth rate with respect to the rotation rate, at least for moderate angular 
momentum. The strong effect of rotation on the growth rate of the spiral modes can thus be 
traced back to this Doppler shifted frequency. What is the mechanism of the instability ? As 
seen in the previous section, the destabilizing role of rotation does not seem related to the 
presence or absence of a coro tation radius, thus discarding a Papaloizou-Pringle mechanism 
(IGoldreich fc Narayanl Il985l ) . Two possibilities have been propos ed for the mechanism of 



SASI without rotat ion; one is the advective-acoustic mechanism ( Foglizzo fc Tagger 200ol: 



Foglizzdl200ll . 120021 ) and the other is the purely acoustic mechanism (IBlondin fc Mezzacappa 
20061 ). Up to now, there is no satisfactory direct argument for the mechanism of the modes 
with a long wavelength. FGSJ07 used a WKB approximation to prove that the instability of 
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the modes with a short wavelength is due to an advective-acoustic mechanism and extrap- 
olated this conclusion to the modes with a long wavelength, which are the most unstable. 
This method, recalled in Appendix B, is based on the identification of acoustic waves and 
advected waves at a radius immediately below the shock surface, and the measurement of 
their coupling coefficients, above this radius due to the shock, and below this radius due to 
the flow gradients. These coupling processes are responsible for the existence of two cycles, 
namely a purely acoustic cycle characterized by an efficiency 71, and an advective-acoustic 
cycle characterized by an efficiency Q. 

By using the same method, the present study does not address directly the instability 
mechanism of long wavelength modes. However, the WKB approximation enables us to de- 
scribe, in a conclusive manner, the instability mechanism of short wavelength modes affected 
by rotation. First we checked that when the shock distance is increased (r s h = 20r*), the 
overtones are also unstable and their growth rate is an oscillatory function of the frequency 
similar to Fig. 7 of FGSJ07. The effect of rotation on the advective-acoustic cycle is il- 
lustrated by Fig 3, for the spiral modes m = ±1 corresponding to the 10-th overtone, as 
a function of the rotation rate. The cycle efficiency Q is strongly amplified by rotation if 
m > 0, while strongly damped if m < 0. The stabilization of the counter- rotating spiral 
coincides with a marginally stable cycle Q ~ 1. The calculation of the amplification factor 1Z 
of perturbations dur ing each purely acoustic cycle indicates its stability (TZ < 1). Contrary 



to the expectation of lLaming I (120071 ) (see next subsection), rotation clearly favours the spiral 



mode of the advective-acoustic cycle. 

This consequence of rotation established unambiguously for short wavelength perturba- 
tions is identical to the influence of rotation on the fundamental mode of SASI: we consider 
this a new hint that the advective-acoustic mechanism can be extrapolated to low frequencies. 
The detailed analysis of the consequences of the Doppler shifte d frequency on the increase o f 



the advective-acoustic efficiency Q will be presented elsewhere (jYamasaki fc Foglizzo I (120081 ). 
in preparation). 



4.2. Comments on the Results of Laming (2007) 



The effect of rotation on the growth rate of SASI, esta blished in Sect. 3 in a cylindrical 
geometry, is qualitatively similar to the effect conjectured by lLaming I (120071 ) (hereafter L07). 
Nevertheless, their investigation about the instability mechanism led them to a different 
interpretation of the roles of the acoustic and advective-acoustic cycles. 



We must point out a fundamental difference between the method of L07 and ours: by 
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using a WKB approximation, we have carefuly defined the range of validity of our method, 
namely short wavelength modes. This guarantees that the advective-acoustic interpretation 
of the instability mechanism is physical and robust, at least in some parameter range. In 
contrast, the existence of a purely acoustic instability is still a conjecture because the domain 
of validity of the method used by L07 is ambiguous: their analytical derivation of a dispersion 
relation when advection is included requires to neglect terms of order (v r /ur) while terms 
of order M. are retained. This approximation is not supported by the results of their Fig. 2, 
which indicates that (f r , s h/^ r sh) is comparable to or larger than -M s h for the modes I = 
and 1 = 1. An accurate description of this acoustic mode, even in a simplified set up, would 
be useful to gain confidence in its possible existence. 

In addition to the question of the validity of the approximations used by L07, we find that 
our results invalidate their reasoning concerning the instability mechanism. They proposed 
that the advective-acoustic mechanism would be essential if r^jr* > 10, whereas a purely 
acoustic unstable process would be dominant for small shock radii, and they argued that 
rotation is a key ingredient to discriminate between the two mechanisms. When rotation 
is included, its effect on SASI has been attributed by L07 to a purely acoustic mechanism, 
despite the results of their Table 3. However, their view that rotation cannot possibly 
enhance the growth of the advective-acoustic cycle is clearly incorrect, at least for the short 
wavelength modes (our Fig. 0J. 



4.3. Consequences of rotation on supernova explosions 

The perturbative study of a simple cylindrical configuration has enabled us to cover 
a large parameter space of shock radii and rotation rates, in order to (i) demonstrate the 
linear selection of non-axisymmetric modes, (ii) establish a correlation between the preferred 
direction of the spiral SASI and the rotation of the collapsing core, (iii) identify the advective- 
acoustic mechanism at work for short wavelength spiral perturbations. 

The fact that rotation favours a s piral mode m = 1, 2 in a cylind rical flow seems directly 



connected to the property observed by lBlondin fc Mezzacappal (120071 ) in their 3D simulations 



including rotation. Tracing back the main influence of rotation to the local Doppler shifted 
frequency u — mQ, we may indeed expect a similar destabilization of the spiral modes with 
a positive value of m, a stabilization of the counter-rotating ones, and a comparatively weak 
influence on the axisymmetric modes. 



Even a moderate amount of angular momentum results in a shortening of the growth 
time of SASI through the destabilization of a non-axisymmetric mode. The promising con- 
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sequences of SASI on both the explosion mechanisms and the pulsar kick could thus be 
considerably modified, since they were established on the basis of axisymmetric num erical 
simulations faurrows et alJbood . 120071 : iMarek fc JankalboO'H : IScheck et allbooi l2006h . Our 
study suggests that the effect of rotation on the linear phase of SASI can be safely neglected 
only for slowly rotating progenitors with a specific angular momentum L <C 2tt ■ 10 14 cm 2 /s. 
Although a fast growth of SASI might be helpful to an early shock revival, the dynamical 
effects of a spiral mode m — 1, and even m = 2, on the possible explosion mechanisms are 
not known yet. 

If the direction of the kick were determined by the geometry of the most unstable / = 1 
SASI mode, our perturbative approach would suggest a kick-spin misalignment. The strength 
of the equatorial kick may be diminished by the domination of a symmetric mode m — 2. 
It is worth noting however that the relationship between the timescale of the most unstable 
SASI mode and the onset of explosion is not straightforward, and should be evaluated by 
future 3D numerical simulations. Our linear approach modestly aims at guiding our intuition 
for the interpretation of these simulations. 



A. Derivation of the Basic Equations 
A.l. Basic Equations 

The basic equations describing the flow are 



dp 
dt 



+ V • (pv) = 0, 



dv 
~dt 



+ w x v + V 



\v\ c 



7 



-VS, 



(Al) 
(A2) 

ot " - P ^ 

Small amplitude perturbations are superimposed onto the above equations. From the trans- 
verse components (9, z) of the Euler equation, the vorticity can be expressed as follows: 



7-1 

dS „ n C 



ik z 1 

V r V 




L 

r 


. Svt 

lUJ — 

V r 


>-+ 


im 
rv r 



c 
7 



lUJ 
V r 



5w z 



From the definition of the vorticity vector, 



-5S -5f-5q 



. im 
0Wr = — ov 7 



ik z 8v e , 



(A4) 
(A5) 

(A6) 
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drSvg 

dr 
d5v z 

dr 



imSv r + rSw z , 
ik z 5v r — Swg. 



(A7) 
(A8) 



Using the two equations (1A4j) and (1A5j) and the definition of vorticity leads to: 

4- - — 1 (rSw r ) = 0. (A9) 
dr v r J 

The radial Euler equation combined with Eq. (jSJ) and ([9]) is 

d5/ T ■"' 



dr 



zu;5tv H + 



55 -<Jg 



(A10) 



Guided by the conservation of 5K in a radial flow (jFogfizzoll200ll ). let us define the quantities 
SKi, 5K 2 as follows: 



(All) 
(A12) 
(A13) 



8K\ = v r rSw z — im — 5S — Sq } , 

7 



iujrbvQ — imSf, 



v r Swg + ik z — SS — Sq 

.7 



-iuj'8v z + zA^ ( Sf 8v< 



(A14) 



where Eq. (1A12I) and Eq. ( 1A14I) are deduced from Eqs. (1A4I) and (1A5I) . The flow quantities 
Svg, Sv z , 5wg, Sw z can be expressed with SK\ and SK 2 using Eqs. ( 1A11HA141) : 



5v z 
5 w e 
Sw z 



m Sf %SK X 



r u 
5f 



to w 



Lk 7 , 



k z — + -,[ SK 2 + -^SK, 



ik z 



5q 



7 



-SS) + 



SK 2 



im 



v r r 



[Sq SS 



7 



SKt 



Using Eq. (JATJ), (JA1| and (jATOl) . we can prove that 



dr v r 



d iu 1 
dr v r 



6K, 



2L 



^Sw r . 



(A15) 
(A16) 

(A17) 

(A18) 

(A19) 
(A20) 
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The perturbations of radial velocity, sound speed and density are related to /, h, SS, Sq as 
follows: 

+ (A21) 
c rc c J 

—Jv e + 5 4- MHh - I 
rc c 

—Jv e + 5 4 - M 2 Sh -SS). (A23) 
rc 1 c z 



Sv r 




1 | 


'Sh 


v r 


1 - 


-M 2 1 


Sc 2 


7 


- 1 | 


(5f 


c 2 


1 - 


- M 2 1 




dp 




1 | 


(6f 


P 


1 - 


- M 2 1 


\c 2 



The continuity of mass flux is 



dr v r p rv 
The differential system satisfied by Sf, Sh, SS, Sq is 



dr v r (l-M 2 ) y u c 2 1 1 ^ c 2 

L 8K X 



v r r 2 1 — M 2 



d5h - ' J (^Sf-MHh-SS+ 5 4 



dr v r (l — M 2 ) \c 2 lu c 



(A22) 



dSh iuj' Sp im . , . „ 
Sv e . (A24) 



(A25) 



6K X ( Lu' , 

™ + T7i > ( A26 ) 



ujr 2 v r \ c 2 (l — .M 2 ^ 

17 = ^ S + Hix)' (A27) 

di « ^.j, + if A V (A28) 



dr t> r \ pv r 

The transverse velocity perturbations Svg, Sv z at the shock are expressed in Appendix A. 2. 
by: 

8v e , S h = — A((v rA - v r)Sh ), (A29) 
Sv Zi8h = ik z A((v rjl -v r>sh ). (A30) 

Together with the boundary conditions (Eqs. (I12H15I) ) established in Appendix A. 2., we 
deduce from the definition of Sw r ,SKi,SK 2 that these three quantities vanish at the shock. 
From the conservation Eqs. (1A9jl . (1A19jl . (1A20I) . we conclude that Sw r , SK\ and SK 2 are 
uniformly zero throughout the flow. The flow quantities Svg, Sv z , Swe expressed in Eqs. (1A15I - 
IA18|) are thus simplified accordingly, and the differential system (1A25HA28]) is transformed 
into the simpler Eqs. 
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A.2. Boundary Conditions 



The Rankine-Hugoniot relation is written as 



Pi S IVi 



Piiyx ~ v s ) ■ — , 
|n| 

JL.X + PlCl 

In] 



7 



(vi 



vi - v s 



n 

|n| 



+ 



7 



Psh(v sh - V s ) • — , 

n 



/'si, i (v sh - v s ) • — 
|n| 

(v sh - v s ) • ti, 
(v sh - v s ) • t 2 , 



+ 



7 



7 



(A31) 

(A32) 

(A33) 
(A34) 

(A35) 



where v s is the velocity vector of the shock surface and n, ti and t 2 are the vector normal 
and tangent to the shock surface which is written at first order as follows, 



n 

ti 
t 2 



I 1 dr s dr s 
r s 88 ' dz 



dr s 



fa" 



Considering small perturbations of the above Eqs, flA31[) - (lA35[) . we obtain 
PsbVr^Shab. + iuj'A((p sh - pi) = AC 



^rsh^Psh + 2p sh tV iSh <5iV iS h + -p S hC s h^C sh + 5 p. 



7 



_sh 

1 



AC 



AC 



dr \ 2 7 — 1 J 1 dr 



2 7-l/ sh j 

and Eqs. (jA29ll - (jA30j> . Using the relat ions in the steady flow, 

d pv r 

— (pV r ) = , 

dr r 



(A36) 
(A37) 
(A38) 



(A39) 



(A40) 



(A41) 



(A42) 
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dr V 2 



d_ 

dr 




(Pvr + P) 



r 



(A43) 



(A44) 



/w r dr ' 



we obtain the boundary conditions f|T2l) - f|T5l) . 

B. WKB Method for the Calculation of the Amplification Coefficients 

In order to interpret the complex eigenfrequency u of a given eigenmode, we compute 
the efficiency of the advective-acoustic and purely acoustic cycles associated to the real 
frequency u r = Re(co>) of this eigenmode. The coefficient lZ(u r ) is defined by the amplifica- 
tion of perturbations after one purely acoustic cycle, initiated at the shock by an acoustic 
wave propagating downward. Q{uj r ) and Q q (u r ) measure the amplification of pertubations 
through an advective-acoustic cycle, initiated at the shock by the advection of an entropy 
perturbation 5S with 5q = 0, or a heat perturbations Sq with SS = respectively. Each of 
these coefficients 71, Q and Q q is the product of the coupling coefficient at the shock {TZ^, 
Qsh, or Q q sh ), multiplied by the coupling coefficient through the flow (TZy, <2v or 0,%). The 
technique of calculation of each factor is the same as that described in the Appendix D of 
FGSJ07. The calculations are based on the decomposition of the variables onto the basis of 
advected and acoustic perturbations, which is exact when the flow is uniform. Even when 
the flow is moderately inhomogeneous, a similar decomposition is obtained using a WKB 
approximation. Since the definitions of the variables 5f, Sh, SS, 5q employed in this paper 
are slightly different from those in FGSJ07, the decomposition is modified as follows: 



The superscripts +, — , S, q refer to the contributions of the ingoing and outgoing acoustic 
wave, the advected quantities SS and 5q, respectively. Adopting the WKB approximation, 
the quantities Sf ± ,8h ± associated with the acoustic waves satisfy the differential system 
(J6H9J) where 5S = and 5q = 0, and the radial derivatives are replaced by a multiplication 
by ik±: 



6f = 5f + + 5f- + 6f b +6P, 
Sh = 5h + + 5h~ + 5h s + 5h q . 



(Bl) 
(B2) 




(B3) 



(B4) 
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The dispersion relation of acoustic waves corresponds to: 

J MTU m ,x 

k± = -Y^W (B5) 

The advected quantities satisfy the differential system (JHHSJ) where the radial derivatives are 
replaced by a multiplication by iko, where k Q = to' /v. 

ikoSf = \mW - M^ 6 -4 + [1 + (7 - l)^ 2 ]- ) , (B6) 

tv(i — M. 2 ) y uj c 2 7 J 

ik Sh s = *^ (£-5f 3 - M 2 5h s - SS) , (B7) 
v r (l — M. 2 ) \c z uj J 



vJl — M 2 ) \ u> c 2 c 2 1 



ikoSh* = lU ' (^-Sf g - M 2 5W + *J ) . ( 09 ) 

v r (l — M. 2 ) \C 2 UJ c 2 
Solving these two sets of equations leads to: 



l f> = HL^-t ls , (Bii) 



uj' 1 — fi 2 M 2 7 

S h s = ^ S f s -5S, (B12) 
uj c 2 

ft. = (B14) 
1 - \i 2 M 2 c 2 

The coupling coefficients TZ s h, Qsh and Qg h are obtained by decomposing the variables at the 
boundary described by Eqs. (I12H15|) onto the basis of acoustic and advected perturbations, 
immediately below the shock: 

SU = Sfc + Sfc + Sfl + Sfa, (B15) 

n sh = (Bie) 

Qsh = (B17) 



S fsh 

s 



Qsh = (B18) 
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The three coupling coefficients TZ\/, <2v and <2y are calculated by measuring numerically, at 
a radius R immediately below the shock (R = r s h), the acoustic feedback 5f~(R) that would 
be produced, either by an ingoing purely acoustic perturbation Sf + , or a purely advective 
perturbation 5f s , or 5f q . Each of these three coefficients is calculated by integrating the 
differential system (JMl]) from the radius R down to the accretor surface. For example, the 
boundary condition used at r = R for the calculation of Qv involves a perturbation of 
entropy and vorticity Sf s (R), and the right amount of acoustic feedback 5f~(R), 

5f(R) = 5f(R) + 5f-(R), (B19) 

such that the inner boundary condition at the accretor surface is satisfied. The coupling 
coefficient Qv measures the efficiency of this acoustic feedback: 

Since Q g is negligible compared to both Q and 71, we discuss only 1Z and Q in the text. The 
WKB decomposition is a good approximation when the inhomogeneity caused by the conver- 
gence of the flow, gravity and cooling is moderate within a wavelength of the perturbation. 
Since we use this decomposition immediately below the shock front, the approximation is 
valid when the inhomogeneity of the steady flow just below the shock is sufficiently small. 
The amplification coefficients Q and 1Z illustrated in our Fig. H] were computed in a flow with 
a large shock radius (r s h = 20r*), for a short wavelength mode (tenth overtone), in order to 
obtain reliable results. 
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L 

Fig. 1. — Growth rate and corotation radius r co /r* as a function of the specific angular 
momentum L, when r sh /r* = 5. The growth rate is normalized by |f r , s h|/( r sh — r *)- L is 
normalized by 27r-10 12 [cm 2 /s], and corresponds to the rotation frequency / p [Hz] extrapolated 
at 10km. Thick lines are for the modes with m > and thin curves are for m < 0. The 
corotation radii are displayed only for the modes with m > 0. The solid, dotted and dashed 
lines are for the fundamental modes, first and second overtones respectively. Left: The case 
with m — ±1. Right: The case with m — ±2. 
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Fig. 2. — Growth rate of the axisymmetric mode (m = 0) as a function of the specific angular 
momentum L, when r sh /r* = 5. The growth rate is normalized by |f r , s h|/( r sh — r *)- L is 
normalized by 27r-10 12 [cm 2 /s], and corresponds to the rotation frequency f p [Hz] extrapolated 
at 10km. The solid, dotted and dashed lines are for the fundamental modes, first and second 
overtones respectively. The value of the wavenumber k z is indicated on each plot. 
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Fig. 3. — Azimuthal wave number m of the mode with the largest growth rate, for each 
specific angular momentum L and shock radius r sh . Units are identical to Fig. 1. In the 
dotted area, the transition between m — 1 and m = 2 is irregular. The spiral mode m — 1 
is unstable in most of the parameter space (above the dashed line). 
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Fig. 4. — In a flow with r s h/r* = 20, the two spiral components m = +1 (thick lines) and 
m = — 1 (thin lines) of the 10-th overtone are analyzed as a function of the specific angular 
momentum L, by solving the boundary value problem (linear scale) and computing the cycle 
efficiencies <2, 1Z in the WKB approximation (logarithmic scale). The growth rates Im(cj) 
(full lines) are compared with the advective-acoustic efficiencies Q (dashed lines, logarithmic 
scale) and the purely acoustic efficiencies 1Z (dotted lines). Units are identical to Fig. 1. 



